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. . Abstract 

Q 

0^ ■ It is well known that the appearance of non-reversibility in classical chaotic systems is connected 

, with a local instability of phase trajectories relatively to a small change of initial conditions and 

parameters of the system. Classical chaotic systems reveal an exponential sensitivity to these 

^T) • changes. This leads to an exponential growth of the initial error with time, and as the result after 
> . 

. the statistical averaging over this error, the dynamics of the system becomes non-reversible. 

(N 

' In spite of this, the question about the origin of non-reversibility in quantum case remains 
OO ! 

' actual. The point is that the classical notion of instability of phase trajectories loses its sense 



o 



> 



during quantum consideration. 

The current work is dedicated to the clarification of the origin of non-reversibility in quantum 
chaotic systems. 



^ ' For this purpose we study a non-stationary dynamics of the chaotic quantum system. By analogy 

H : 

with classical chaos, we consider an influence of a small unavoidable error of the parameter of the 
system on the non-reversibility of the dynamics. 

It is shown in the paper that due to the peculiarity of chaotic quantum systems, the statistical 
averaging over the small unavoidable error leads to the non-reversible transition from the pure 
state into the mixed one. 

The second part of the paper is dedicated to the kinematic description of the chaotic quantum- 
mechanical system. Using the formalism of superoperators, a muster kinematic equation for chaotic 
quantum system was obtained from Liouville equation under a strict mathematical consideration. 
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INTRODUCTION 



Study of the quantum reversibility and motional stability is of great interest [1-4] . This 
interest is due to not only the fundamental problem of irreversibility in quantum dynamics, 
but also to particular application. In particular, it reveals itself in relarion to the field of 
quantum computation. A quantity of central importance which has been on the focus of 
many studies [5-9] is the so-called fidelity f{t), which measures the accuracy to which a 
quantum state can be recovered by inverting, at time t, the dynamics with a perturbed 
Hamiltonian 

m^Me^"'e-^""'m'. (1) 

Here ip is the initial state which evolves in time t with the Hamiltonian Ho, while H = 
Hg + V is the perturbed Hamiltonian. The analysis of this quantity has shown that under 
some restrictions, the series taken from f{t) is exponential with a rate given by the classical 
Lyapunov exponent [5-6]. But here a question appears. The point is that the origin of the 
dynamic stochasticity, which is the reason of irreversibility in classical case, is directly related 
to the nonlinearity of equation of motion. For classical chaotic system this nonlinearity leads 
to the repulsion of phase trajectories at a sufficiently quick rate [10-13]. In case of quantum 
consideration, the dynamics of a system is described by a wave function that obeys a linear 
equation and the notion of a trajectory is not used at all. Hence, at first sight it seems 
problematic to find out the quantum properties of systems, whose classical consideration 
reveals their dynamic stochasticity. For this reason, according to the widely accepted opinion 
the notion of "Quantum Chaos" includes phenomena, related to the quantum- mechanical 
description of systems chaotic in the classical limit [14-17]. 

Unfortunately, our understanding of quantum dynamics of chaotic systems is still quite 
limited. The majority of the existing quantum chaos literature concentrates on understand- 
ing the properties of eigenf unctions and eigenvalues. 

In this paper we consider the issue of irreversibility for quantized chaotic system. For 
simplicity we shall discuss a periodic driving. 

In given work we shall try to show that in essentially chaotic domain, quantum chaotic 
dynamics is characterized by the transition from the pure quantum-mechanical state into 
the mixed one. In this case, irreversibility in the system appears as a result of loss of 
informaintion about the phase factor of the wave function. To prove this phenomenon, we 
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shall consider model Hamiltonian studied in [1,18]. We shall try to give our explanation 
of the mixed state formation in a chaotic quantum-mechanical system. The main purpose 
of given work is the deriving of kinetic equation for chaotic quantum-mechanical system. 
Interesting method of how to derive a kinetic equation for semi-classical chaotic system was 
offered in [16]. In this work we shall try to do the same in exceptionally quantum domain, 
without application to the semi-classical methods. The paper is organized as follows. In 
the first part formation of a mixed state is considered. The second part is dedicated to the 
kinetic description of quantum chaos. 

FORMATION OF MIXED STATE. 

During the discussion of nonreversible evolution of the quantum- mechanical system nat- 
urally the question emerges. How does a non-reversibihty originate in a quantum system? 
If the quantum system evolves according to the Schrodinger equation, how can a pure 
quantum-mechanical state become a mixed one? The question is that, in contrast to the 
classical chaos, quantum-mechanically irregular motion cannot be characterized by extreme 
sensitivity to tiny changes of the initial data. Due to the unitarity of the quantum dynamics, 
the overlap of two wave functions remains time-independent | (©(i) P = | (0(O)|.^(O)) p, 
provided time-dependence of Q(t) and ^(t) is generated by the same Hamiltonian. However, 
an alternative characterization of classical chaos, extreme sensitivity to slight changes of the 
dynamics does carry over into quantum mechanics. 

As it was mentioned in the introduction, the model Hamiltonian, studied in [1,18] is taken 
as a model for our considerations 

H{Q, P, x) = ^(P2 + p| + QI + Ql) + :,Q\Ql (2) 

where Pi, P2, and Qi, Q2 are canonical momenta and coordinates. Here and further we 
shall consider dimensionless quantities. 

It follows from works [1,18] that for values of parameters x = ^ 1 and £' fsi 3 (where 
E is energy of system), dynamics related to Hamiltonian (2) is chaotic and is characterized 
by correlation function, the width of which is equal to Tc = 1 [1] . Our goal is to study mech- 
anisms of irreversibility and mixed state formation in non-stationary quantum dynamics. 
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For this let present Hamiltonian (2) in the following form 

H{t)^Ho + V{t), 

Ho = ^{P^ + P| + Ql + Ql) + XoQlQl, (3) 

Vit) = VoiQ)foit), Vo = Ax„ • QlQl 

As one can see from (3), the time dependence of the Hamiltonian H{t) is provide by 
amplitude modulation of the parameter x — Xo + Axo/o(t), where fo{t) is a periodic function 
of time with period. Hereinafter we shall work in a domain of chaotic motion Xo ~ 
1, Axo < Xo- 

After taking (3) into account, the Schrodinger equation for the wave function takes the 
form 

.dmMm ^ H{Q,^^)mQ^xm- (4) 

The solution of the time-dependent Schrodinger equation (4) can be written formally 
with the help a time-dependent exponential [14] 

U{t) = exp[-i f dt'H{t%, 

Jo 

where the positive time ordering requires 



[A(t)B{t% 



A{t)B{t') if t>t' 
B{t')A{t) if t<t' 



In our case H{t + Tgh) = H{t), k = 1, 2..., the evolution operator referring to one period 
To , the so-called Flouqet operator U{To) = F [14], is worthy of consideration, since it yields 
a stroboscope view of the dynamics 

mkTo)) = (Prim)- 

The Flouqet operator being unitary has unimodular eigenvalues. Suppose we can find 
eigenvectors \(p^) of the Flouqet operator 

M'Pe) = <^xe- (5) 
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Then, with the eigenvalue problem solved, the stroboscopic dynamics may be written ex- 
phcitly [14] 

lV'n(A:ro)) = Ee-^'^^(^xl^n(0))|9^,). (6) 

X 

Here |'0n(O)) is the eigenfunction of the Hamiltonian H (3), corresponding to the value 
of parameter x(0) — Xq-\- Axo/(0). 

As it was mentioned above, our aim is to prove that one of the signs of the emergence 
of quantum chaos is a formation of the mixed state. Being initially in a pure quantum- 
mechanical state, described by the wave function \ipn)i the system during the evolution 
makes an irreversible transition to the mixed state. 

The information about whether the system is in the mixed state or in the pure one, 
may be obtained from the form of the density matrix [19-21]. Introducing the definition 
Cnm{k) — e~^'^'^"' {(f)jn\'il^n{^)) : let US rewrite (6) in more convenient way for further usage: 

\MkTo))^J2Cnm{k)\<f>m). (7) 
m 

Next, according to the commonly accepted rules [22], the density matrix of the system may 
be expressed via the expansion coeficients Cnm{k) 

Pnm{k)^Y.Cnp{k)Cl^{k). (8) 
P 

After the substitution of the explicit form of Cnm{k) coefficients into (8), we get: 

Pnm{k^ ^mn^- ^'^ \ (9) 

where 

Amn = T.{<l>m\M^)){M^)\<l>n)- (10) 
P 

As seen from (9), non-diagonal matrix elements of the density matrix contain fast oscil- 
lating in time exponential phase factors e^^^'^rn-<t>n) qj^^]^ ^j^g diagonal matrix elements 

Pnn = ^nn = EK'/'n|^p(0))|^ (11) 
P 

Exponential phase factors of the non-diagonal matrix elements Pnm{k) express the prin- 
ciple of quantum coherence [23] and correspond to the complete quantum-mechanical de- 
scription of the system in pure quantum-mechanical state. While they are not equal to zero. 
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the system is in the pure state, and zeroing of these elements is the sign of transition into 
mixed state [20-22]. 

But here a question appears: if the system was initially in the pure quantum-mechanical 
state, how can the transition into mixed state take place? Before answering on this question, 
let us recall, that existence of non-complete quantum information about the state of the 
system corresponds to the mixed quantum-mechanical state [19-23]. I.e. for transition 
into the mixed state, a partial loss of information about the state of the system must occur. 
According to the commonly accepted notions [20] , this corresponds to the loss of information 
about the phases of the system and to the zeroing of non-diagonal elements of the density 
matrix. 

So, to prove the formation of the mixed state one has to show the zeroing of non-diagonal 
elements of density matrix, the equality of which to zero is a sign of a mixed state [19-21]. 

Let us recollect that value 0„ is the eigenvalue of the Floquet operator. Owing to the 
non-integrability of the system, eigenvalues can be obtained only by way of numerical di- 
agonalization of the Hamiltonian Hq (3). According to the main hypothesis of the random 
matrix theory [14,17], the elements of this matrix in the chaotic domain are random numbers. 
So, their eigenvalues can also be considered as random values. 

This statement is valid only in the chaotic domain [24,25]. According to this, the phase 



in exponential factors of the non-diagonal matrix elements in (9) is random quantity with 
all its properties. Such an interpretation of the eigenvalues of chaotic quantum-mechanical 
system is stated in recently pubhshed works [24,25]. In these papers the analogy between 
spectral characteristics of a Hamilton system and random time series is discussed. Further 
we shall treat f{n,m) phase as a mathematical random variable putting aside its physical 
meaning. According to the preceding, it is clear that the values of matrix elements of density 
matrix of the chaotic quantum- mechanical system pnm{k) are random values too. Taking a 
statistical average of expression (9), we have: 



where (...) denotes the averaging over the statistical ensemble. 

Before proceeding the discussion, we should define what one implies under quantum- 
statistical ensemble in our case more precisely. 



f{n,m) = 0^-0, 



(12) 



(13) 
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Recall, that while studying classical chaos, one usually examines the stabihty of the 
system relatively to a small change of initial conditions and system parameters. A small 
initial error of these parameters always exists and remains unavoidable (one may measure 
the parameters of the system at a very high precision, but even in this case there is still a 
small error, the removal of which, i.e. the measurement at absolute precision, is impossible 
[10]). So that, not the existence of the unavoidable error is fundamental, but what kind 
of influence it brings over the system dynamics. It is well known that in case of regular 
systems such an influence is negligible, but if the system is chaotic, the effect of it increases 
exponentially [12]. 

Making a complete analogy to classical chaos, let us assume, that the parameter, charac- 
terizing interaction between oscillators Xq, has some small unavoidable error 6xq. Then the 
set of chaotic Hamiltonians (3), with the values of interaction parameter taken from interval 
[xq — Sxq, Xq + Sxq], may be considered as quantum-mechanical ensemble. According to the 
random matrix theory [14,17,27], there is a random distribution of distances between levels 
that corresponds to ensemble of chaotic Hamiltonians. In particular for chaotic Hamiltoni- 
ans with real matrix elements the distribution of the distances between levels obeys normal 
Gaussian distribution [14,17]. 

The further study of (13) can be performed using more rigorous mathematical substan- 
tiation. For this let us recall some details from the probability theory [26] . 

a) In general case, under the characteristic function of the random variable X, mathe- 
matical expectation of the following exponent is meant 

F{t) = M(exp(itX)), (14) 

where t is a real parameter. 

b) Mathematical expectation itself is defined as a first initial moment /ii of the random 
variable X 

M{X) = {X)=fi, = J2xkPk, (15) 

k 

where X is the discrete random variable which takes possible values xi,X2, ■■■ with appro- 
priate probabilities Pk — P{X — Xk), < ■■■ > means average. 

Taking (14), (15) into account and considering f{n,m) as a random value, we get: 

{Pnm{k))^A^^F{k), (16) 
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where 

F{k) = (e^'=/("'"^)), (17) 

is the mathematical expectation of the characteristic function of the random phase /(n, m), 
and 

Amn^{(f>m\^m{ll^m<Pn). (18) 

Let us assume that random phase f{n,m) has a normal dispersion [26]. This assumption 
is based on the fact that matrix elements of Hamiltonian (3) are real values. Then from (17) 
one can obtain 

F(A;) = e^«V^, (19) 

where 

a^M{f{n,m)) (20) 
is the mathematical expectation of the random phase, 

a'^M(f(n,m))-(M(f(n,m))r. (21) 

On basis of obtained expressions (16), (19) one may say the following. For small k, the 
system is still in a pure state (non-diagonal matrix elements differ from zero). But after 
some time the system passes to the mixed state as the non-diagonal matrix elements of the 
density matrix decrease with exponential law 

Pnm(A:)~e-^. (22) 

This phenomenon is connected with the "phase incursion" [12]. Uncertainty of phase in (9) 
is accumulated little by little with time. Finally, when k > V^/cr, the uncertainty of phase 
is of order 27r, and the phase is totally chaotized (see [3,4]). As the result the system passes 
into mixed state. 

But here a question emerges. How does the uniqueness of the nonlinear chaotic system 
show itself? If we assume the existence of initial dispersion of parameter for linear sys- 
tem, shall we obtain the same result? Before answering on this question let us consider a 
pecuharity of already obtained result (22): 

When approaching to (22), the fact of randomicity of the distribution of distances between 
levels of energy terms was fundamental. We should notice, that this randomicity was made 
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conditional on the randomicity of the distribution of energy termes themselves and not due 
to the widening of the hne 50„ of the energy terms 0„ [14,17]. 

Now let us consider a linear system and see whether this condition is carried out or not. 

Consider ensemble of linear oscillators 

that are characterized by a small dispersion of parameter 5(3 <^ {3. 

The spectrum of linear oscillator is well known and has a regular character [20] 

K(a;)=a;(n+^), (24) 



where uj — P/m is the distance between the levels of the system. 

Easy to see, that if uncertainty of parameter /3 is small 5/3 <^ /?, then this may lead to 
the widening of the levels, but not to the mixing of them. As a result, the value of such 
widening in uj units is: 

5En^5u{n+]^), (25) 

where 5a; ^ ^ = UjS. 

It is already clear, that the condition is not fulfilled. In this case the damping of non- 
diagonal matrix elements of density matrix may be caused only by the widening of energy 
terms and not by their random distribution. One should take into account that the widening 
of term depends linearly on the dispersion 5(3, and each of the Hamiltonians makes an equal 
contribution to the average value. This happens because in the case of linear oscillator, 
each of the Hamiltonians of the ensemble has the same spectrum as others have, and the 
averaging is performed over the widening of energy terms. In case of chaotic ensemble, each 
of the Hamiltonians of the system has its own spectrum (due to exponential repulsion of 
the energy terms of the chaotic ensemble [14,17] they are exponentially sensitive to a small 
dispersion) and the averaging is performed over these spectrum characteristics, and not over 
the widening of separate line. When taking such average, the contribution from separate 
Hamiltonians of the chaotic ensemble is not the same and is described by Gaussian esemble 
[14,17]. This is the fundamental and not the formal difference between these two cases. 

Taking all this into account, we obtain for ensemble average matrix elements of linear 
oscillator: 

1 r^i^ 

{pLm = ^ J_^^ pLit, Sp')d{Sp'). (26) 



Considering (24), (25) and (26) we get in the end: 



5/3t\n — m 



2 



,i{m—n)uit 



COS 




— {n — m)ujt 



(27) 



Let US start the analysis of the obtained result. Firstly it is clearly seen, that time 
dependence is not of Gaussian type, but is inverse proportional function of time 1/t. This 
as minimum decreases the rate of formation of mixed state in comparison to the case of 
chaotical ensemble (see Eq.(22)). More than this, the time of formation of the mixed state 
is completely defined by the value of dispersion 5/3: 



Actually this means the following. Observing the evolution of the system on an arbitrary 
large interval of time t e [0, T], one can always select such an accuracy while measuring the 
system parameter 5/3o to fulfill the condition T < Tc. This means that the system will be 
in a pure state for given observation interval t E [0,T] and given accuracy of measurement 



Hence in some sense we have an analogy with classical chaos. A small initial dispersion 
of values has an exponential effect, which reveals itself in formation of mixed state. In case 
of regular system, one may neglect this dispersion with an accurace up to the widening of 
the energy terms. 

Finally let us discuss one more question. Can one formally obtain the result (22) in a case 
of linear oscillator? From a formal mathematical point of view it is possible, but this case 
will not have a physical sense. For this one has to choose an ensemble of linear oscillators 
with predetermined random frequences that obey Gaussian distribution. Next one does the 
averaging over this ensemble, taken over the random frequences and not over the widening 
of lines. But this case will not have anything common with a real physical situation. The 
thing is that the Gaussian ensemble of Hamiltonians in case of chaotic system as well as the 
ensemble of harmonic oscillators in a case of regular system is directly connected with an 
unavoidable error of measurement of physical parameters of the system, and the averaging 
over them is not a pure formality. 

Zeroing of the non-diagonal elements of density matrix is the sign of quantum chaos 
beginnings and formation of mixed state. From that moment the quantum dynamics is 
non- reversible, since the information about the wave functions phase is lost. 
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FIG. 1: The graph of the dependence pnm (9) on time t = kTo . The graph is obtained using 
(9) and by way of numerical diagonalization of the Hamiltonian (3) H{P,Q,x{0)). As it is easily 
seen from the graph that when t = kTg = \/2/a > Tc, Tg = ^ = O.Itc, Tc = 1 [1], zeroing of 
Pnm happens. In the numerical calculations we have used definite initial state ■0(0), where ■0(0) 
is the eigenvector of the Hamiltonian H{P,Q,x{0)),x{o) = Xg + Axo, and definite eigenfunctions 
\4>n), \4>m), n = 6, m = 80 of the Hamiltonian H{P,Q,xo). The dimension of the diagonalized 
matrix was 1000 x 1000, Axo = 0.1. For the statistical averaging of the numerical data we have 
used ensemble of 50 Hamiltonians, corresponding to the small dispersion of parameter Sxq = 0.01. 



After the formation of a mixed state the quantum-mechanical description loses its sense 
and we have to use a quantum-statistical interpretation. Along with this, the notion of 
fidelity loses its meaning, i.e. by the time of the beginning of reverse evolution the system is 
no longer in the pure state described by the wave-function. The reverse transition from the 
mixed state to the pure one and reconstruction of the wave-function are impossible because 
of the irreversibility of the process. 

Given above reasoning can be checked by numerical diagonalization of the Hamiltonian 
(2), defining eigenvalues (p^ and by estimation of the averaged non-diagonal elements of the 
density matrix according to (9). 

The result of numerical calculations is represented on Fig.l. 
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KINETIC DESCRIPTION OF THE CHAOTIC QUANTUM-MECHANICAL SYS- 
TEM. 



In previous section we studied mechanism of mixed state formation. After formation of 
mixed state quantum-mechanical consideration loses meaning and there is a need to use a 
kinetic description. Kinetic equation for the chaotic quantum-mechanical system for the first 
time was obtained in [16]. But this study was done in the semi-classical domain. Namely for 
example, zeroing of non-diagonal part of density matrix was proved by use of semi-classical 
approximation. The main purpose of given work is the deriving of kinetic equation for 
chaotic quantum-mechanical system in exceptionally quantum domain without application 
to the semi-classical methods. For obtaining of muster equation for the density matrix we 
shall use a method of projection operator [28]. 

Let us split the density matrix operator p on a slow pR and a fast varying pisiR operators: 
P — Pr + Pnr- Relevant part pR in the basis \(pn) contains only diagonal elements, whereas 
non-relevant part pnr contains only non-diagonal elements. These elements, as it was shown 
in previuos section, contain fast oscillating exponents and when taking average over the 
ensemble, the zeroing of them takes place. Elimination of the diagonal part from the density 
matrix is a linear operation, which satisfies the property of projection operator — D [28] 



Let us note that this reflection is non-reversible. Due to the zeroing of non-diagonal part 
of the density matrix, part of information is lost. 

Inasmuch as relevant statistical operator pnit) is different from the total operator 
generally speaking it does not suit the Liouville-Fon Neumann equation [21] 



PR = bp, Pnr = (1 - b)p. 



(29) 




dpNR 



+ i{l-b)L{pR + pNR)^0. 



(31) 



dt 
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For the purpose to obtain closed equation for pR we exclude from the first equation (31) 
Pjvr . As a result we get 

+ iDLpRit) + jK{t- ti)pR{ti)dti = -iDLexp[-i{t - to){l - b)L]pNR{to), 

to 

K{t - h) = DLexp[-i{t - ti)(l - D)L]{1 - D)L. (32) 

This equation is valid for t — V^/a > Tc- Evidently pR{t) is expressed by way of values of 
PR{ti) taken for the time interval tg < ti < t, and additionally trough the value of PNRito)- 
If in the initial moment of time t — to the system is in a pure quantum-mechanical state 
then pNRiO) ^ 0. 

In this case solving the equation (32) is problematically. However let us recollect that for 
to > V^/cr system is already in a mixed state. Therefore equation (32) takes more simple 
form {pNR{to) = 0) 



dpRjt) 
dt 



+ iDLpR{t) + 1 K{t- ti)pR{ti)dti = 0. (33) 



For solving equation (33) we shall use a method of super-operators [29]. But before 
performing this, we should notice, that we have obtained a closed equation for relevant part 
of statistical operator. Wc were able to come to this only because when t > to = V^/a the 
system is in the mixed state and all non-diagonal matrix elements of the density matrix are 
equal to zero. 

Further when studying the evolution of the system we shall consider as the origin of time 
the moment of the formation of the mixed state in the system. This corresponds to a formal 
transition to hmit to — > — oo. Further for simphfication of (33) we shall use Abel theorem 

[26] 

hm i r f{t)dt = /(O) - hrn /° e<''^ ^f{t')dt' . (34) 

T->oo I Jt e^+Q J-oo dt' 

Taking (34) into account, Eq.(33) will have the following form: 

dpR{t) 



dt 



+ iDLpR(t) = lim / e<*'-*^K{t - t')pR{t')dt' . (35) 

*+0 J —oo 



According to the method of superoperators [29], the correspondence of one operator to 
another may be considered as representation. The operator in this case will be represented 
by a matrix element with two indices, while the linear product of operators is a matrix with 
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four indices, i.e. a superoperator [29]. A concrete example is the projection of D operator 
on diagonal elements. We should notice, that in our case, the projection operator is some 
definite physical procedure of averaging matrix elements of the density matrix over Gaussian 
chaotic ensemble. 

From the relation Dpnmit) = Pnm{t)Smn we come to the following representation of D 

superoperator Dnmn'm' = Snn'Smm'^nm, SO that 

^ ^ Dmnm'n'Pm'n' Pnm^nm- (^^) 
m'n' 

Taking (36) into account, Eq.(35) is 

+ ib{Lpn)nn = - dt'[K{-t')pR{t + t')Ue'r (37) 

In expression (37) and further for short we shall omit R index for diagonal matrix elements 
of the operator p. Considering the relation 

[L, PR]nn — X/ ^nnaaPaa — 0, (38) 
a 

and representing Liouville operator as L = Lq + L' in compliance with (3) we get: 



dpnnit) 



f dhe<'-'^^Y.^nnmm{t-h)Pmn{tl), (39) 



dt 
and 

Knnmmit) = [L'e-^*(^-^)^(l - L')L']_. (40) 

From (40), (3) and from representation of Liouville operator in form L = Lo + L', one can 
see, that kernel K^nmrniS) is at least of second order by Axq. Next it is easy to check the 
corectness of the expression: 

L'abmm. — T^XKim^bm — Knb^am) = Ki6 — ^ab — 0, (41) 
m m 

for Liouville superoperator 

The relation (41) in its turn leads to the rule of sums Knnmm — 0. Taking the symmetries 
Labcd = Lcdah, Dabcd = D^dab into accouut, we get from (39): 

dpnnit) ^_ ^^^ge(*-ii) ^ [K^^^^{t - h)pmm{tl) - K^mnnit - h)Pnnitl)]. (43) 
J OO i_ 
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Next we shall make the following approximations. With accuracy up to the value of (Axq)^ 
order in the exponent in expression (40) we shall replace the complete Liouville operator 
L — Lq + L' with Lq. With the same precision we may set Pnn{t — t') = Pnn{t). As a result, 
from (43) we arrive at 

dpnn{t) 



J2[^nmPmm(t) -WmnPnn(t)], (44) 



where 



W 



- f dt'e^*'[L'e^*(^-^)^°(l - i))L']_. (45) 
Since 

{DLo)abcd = Sab[(l)aSacSbd " (f^bSbd^ac] = 0, (46) 

in expression (45) the D operator in the argument of exponential function may be omitted. 
Taking into account time dependence of the operator x{t) = Xg + Axof{t), f{t) — 

oo . 

J2 , fl — 27: /To for non-diagonal matrix elements we get 

1^=— oo 



= E (W^mPmmit) - W^^„Pnn(i)), (47) 

oo 

where — f | V^mP J2 S(Enm — i^^) is the transition amplitude between the eigenstates 
of the Hamiltonian Ho (3), Enm = <Pn — 0m , Vnm is the matrix element of the operator 
Vo = AxoQlQl i'^ the basis of eigenfunctions of the Hamiltonian Ho, Vnm = ('^nlKlV'm)- 

Equation (47) describes non-reversible evolution of the system from non-stationary state 
to the stationary state defined by the principle of detail equilibrium [20]. To prove irre- 
versibility of the process let us consider time dependence of non-equilibrium entropy [30] 

S{t) = -KBY.Pnn{t)Hpnn{t)), (48) 
n 

where Kb is the Boltzmann constant. Taking into account J2n Pnn{t) — 1 , from (48) we get 
dS{t) 



dt 



= -Kb E E ^nm[pmm{t) - Pnn{t)] ln(p„„(t))- 



-i^BE%^ = lKBY.Y.Wnrn[Pnn{t)-PrnrnmHPnn{t))-HPrnrn{m. (49) 

Due to the property of logarithmic function 

{Pnn{t) - Pmm{t)){MPnn{t)) - ln{pmm{t))) > 

15 



we see that ^ > 0. This testifies the growth of entropy during the evolution process. 

More exact estimation of entropy growth may be obtained from the principle of de- 
tail equilibrium [20]. It is evident from the expression of transition probability Wnm = 

oo 

f iKirrtP J2 S{Enm ~ '^^) that Only resonant states E^m = are included in the process. 



Transitions between states lead to the redistribution of the initial probabilities given by 
pnn{to)- As a result after the time being we have stationary distribution of the probabilities 
defined by the principle of detail equilibrium. 

According to this principle number of transitions from the state m into the state n is 
equal to the quantity of reverse transitions 



Taking into account that in our case Wnm — Wmn , for the entropy growth we get 



where the number of levels N, completely included into the process is defined by the spectral 
characteristics of the Hamiltonian (3) and perturbation (according to (47) ). 

The authors express their gratitude to Professor A.Ugulava for valuable suggestions and 
useful discussions. 




nn- 



(50) 




(51) 
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